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Abstract 

The notion that Stimulated Brillouin Scattering (SBS) is primarily 
defined by bulk material properties has been overturned by recent work 
on nanoscale waveguides. It is now understood that boundary forces of 
radiation pressure and electrostriction appearing in such highly confined 
waveguides can make a significant contribution to the Brillouin gain. Here, 
this concept is extended to show that gain enhancement does not require 
nanoscale or subwavelength features, but generally appears where optical 
and acoustic fields are simultaneously confined near a free surface or ma¬ 
terial interface. This situation routinely occurs in whispering gallery res¬ 
onators (WGRs), making gain enhancements much more accessible than 
previously thought. To illustrate this concept, the first full-vectorial an¬ 
alytic model for SBS in WGRs is developed, including optical boundary 
forces, and the SBS gain in common silica WGR geometries is computa¬ 
tionally evaluated. These results predict that gains 10^ times greater than 
the predictions of scalar theory may appear in WGRs even in the 100 pm 
size range. Further, trapezoidal cross-section microdisks can exhibit very 
large SBS gains approaching 10^ m"^W“^. With resonant amplification in¬ 
cluded, extreme gains on the order of 10^^ m“^W"^ may be realized, which 
is 10® times greater than the highest predicted gains in linear waveguide 
systems. 


1 Introduction 

Stimulated Brillouin Scattering (SBS) is a fundamental acousto-optical nonlin¬ 
earity [1-3] that occurs in all states of matter and is a powerful tool for coupling 
traveling photons to traveling acoustic phonons. The process occurs through 
light scattering from photoelastic perturbations created by an acoustic wave in 


1 



a medium, accompanied by simultaneous amplification (for Stokes scattering) or 
attenuation (for anti-Stokes scattering) of the acoustic wave through electrostric- 
tive feedback. When sufficient optical power is pumped into the medium, the 
Stokes scattering occurs regeneratively and a Brillouin laser is generated at a 
fixed offset from the pump laser. In addition to being a method for generat¬ 
ing [4], suppressing [5], and interrogating [6-8] acoustic modes, SBS has been 
employed effectively in a wide range of applications [9] including ultra-narrow 
linewidth lasers [1], quantum optics and electrodynamics experiments [10-13], 
distributed temperature and pressure sensing [14], optical amplifiers [15] and 
modulators, photonic oscillators [16-18], optical cooling [5], optical phase con¬ 
jugation [19], slow light [20,21], and nonreciprocal optical devices [20,22,23]. 
SBS is one of the highest gain material-level nonlinear it ies in optics [15] making 
it one of the first nonlinear effects to appear in high energy density systems such 
as high power fiber lasers and high quality factor resonators. The suppression of 
SBS is thus particularly important in laser systems as it represents a parasitic 
loss effect and prevents further power scaling of the source [24,25]. 

We consider in this study the SBS gain, which describes the amplification or 
attenuation experienced by a small signal probe that is offset from the pump by 
the acoustic frequency. Textbook descriptions [15,26] state that the line-center 
SBS gain is defined by the electrostrictive constant of the material, the optical 
pump frequency, the refractive index, the speed of sound, and the phonon losses 
in the material. The role of geometry or boundaries is generally not consid¬ 
ered. However, in recent theoretical [27] and experimental [28, 29] studies on 
devices with nanoscale dimensions, significant SBS gain enhancement has been 
deduced. This surprising result occurs due to very large radiation pressure and 
electrostriction forces at the device boundaries [27,28,30,31]. Here we propose 
that these boundary force enhancements will not only appear in sub-wavelength 
devices, but also whenever an optical field is confined near a material interface 
or free surface. Such surface-confined Brillouin scattering interactions are al¬ 
ready known to occur in whispering gallery resonators (WGRs) as demonstrated 
in [4,5,16,17,20,32,33]. In microscale WGRs, light and traveling acoustic waves 
are both tightly confined near the resonator surface, even without the benefit 
of any sub-wavelength dimensions, thus very large SBS gains should be ex¬ 
pected. A similar insight has also been developed recently [34] in a study where 
the effects of boundary forces on axisymmetric stationary vibrational modes in 
WGRs were analyzed. Experimentally, the highest measured SBS gain of 4x 10^ 
m'^W'^ has been observed in a microsphere WGR in a Forward SBS configura¬ 
tion (see Supplement Table S.l in ref. [20]), and it is important to understand 
how such large values of gain arise in such large systems. 

In this paper, we demonstrate that the large optical forces generated at 
WGR boundaries can both enhance and attenuate the Brillouin gain by mul¬ 
tiple orders-of-magnitude. We provide, for the first time, a full vectorial de¬ 
scription of SBS gain in WGRs incorporating the effects of boundary forces 
including electrostriction and radiation pressure. We then use the formalism 
to computationally evaluate the SBS gain in silica WGRs of several commonly 
used geometries, across a range of size scales, while examining the individual 
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contributions of the optical forces to the SBS gain. 


2 Theory of SBS in Linear Systems 

SBS is a three-wave interaction between two optical fields and one traveling 
acoustic field in a dielectric medium [2,15]. Coupling between the waves is 
enabled primarily through photoelastic scattering and electrostriction. The 
two optical fields are selected such that they generate a spatiotemporal over¬ 
lap matching the frequency and momentum associated with the acoustic wave, 
resulting in the phase matching requirement described below. This field over¬ 
lap generates a spatially varying distribution of electrostrictive pressure. The 
gradient of this pressure generates optically-induced forces which can enhance 
(or even attenuate [5]) the acoustic wave. The acoustic wave, in turn, creates a 
periodic refractive index perturbation in the medium through the photoelastic 
effect, through which light scattering between the optical fields can occur. The 
higher intensity, typically higher frequency optical field is called the “pump”, 
while the scattered or probed optical field is called “Stokes” due to its lower 
frequency. SBS can also occur using optical signals of two different polariza¬ 
tions [22,35]. When light scattering takes place in the forward direction, MHz 
- GHz frequency acoustic waves are involved, and this is called Forward SBS 
(F-SBS). Scattering in the backward direction is called Backward SBS (B-SBS), 
in which multi-GHz acoustic waves mediate the optical coupling. 

A key requirement in SBS is the energy and momentum phase matching. The 
mathematical relationship for Stokes light scattering is written diS kp — ks = K 
and ujp — ujg = where kp^ kg^ and K are the momentum vectors of the 
photon modes and the phonon mode respectively, while and Vi are the 

corresponding angular frequencies. 

An additional requirement of SBS is that the optical forces be able to effec¬ 
tively drive the acoustic mode; this is analogous to having strong coupling of 
power from the pump beam into the Stokes and acoustic modes. The deriva¬ 
tion in Section 4 is centered on calculating this coupling using the opto-acoustic 
overlap integral (see Eq. (2)) of the form {F\U). Good coupling occurs when 
the optical forces F align well with the acoustic mode [/, which is intuitively 
understood as efficient power transfer into the acoustic mode. An alternative, 
but entirely equivalent, method for determining the strength of this coupling is 
to use the optical coupling integral of the form (Ei | Ae|E 2 ). The optical coupling 
integral describes how a perturbation Ae to the dielectric permittivity (induced 
by the acoustic mode for SBS) facilitates power exchange between two opti¬ 
cal modes Ei and E 2 . The equivalence of these two formulations was recently 
proved in the case of low loss and reversible optical forces (valid for the purposes 
of this paper) by Wolff et al. [36]. Because of these overlap integrals, SBS can 
be symmetry-forbidden for certain optical and acoustic mode combinations (not 
considered here). More detailed discussions of the permitted transitions can be 
found in [30,37]. 

Even though SBS coupling is typically considered as being dominated by 
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bulk optical electrostriction, recent work has demonstrated that boundary op¬ 
tical forces play a significant role in confined geometries. These effects become 
extremely important in nanoscale waveguides [27-29,38]. The electrostrictive 
boundary forces arise because of the difference in optically induced stress across 
the waveguide-cladding boundary, while radiation pressure occurs due to mo¬ 
mentum exchange from photons reflecting off these surfaces. Radiation pressure 
can also manifest within the bulk of the material [39], but this effect is orders- 
of-magnitude smaller than the effects considered here. Better overlap between 
all these optical forces and the acoustic wave in tightly confined waveguides 
increases the SBS gain as well. Recent experimental work [28] on linear waveg¬ 
uides with sub-wavelength cross-sections has demonstrated enhancement of SBS 
gain due to the appearance of these large boundary optical forces and the high 
degree of modal overlap. Based on these insights, mathematical models for 
computing SBS gain in linear systems have already been developed offering ro¬ 
bust predictions for straight waveguides [27,30,36], slot waveguides [38], and 
fibers [35]. The action of these forces in WGRs has remained unexplored. 

3 The Difference Between Waveguides and 
Whispering-Gallery Resonators 

The significant field enhancement within high-finesse WGRs enables very low 
pump power threshold for almost any nonlinear process, from the laboratory 
reference frame. This is because the optical power circulating within WGRs 
is amplified with respect to power coupled in from the waveguide by a factor 
of Finesse/27r. This amplification factor can easily exceed 10^ for microscale 
silica spheres [40], shells [41], and disks [18]. For SBS in WGRs, both the 
pump and Stokes beams are resonant with comparable finesse, so the resonant 
enhancement of gain can be on the order of 10^^ for ultra-high-Q devices. The 
pump finesse contributes due to resonant enhancement of the pump power, while 
the Stokes finesse accounts for the fact that the optical path length is finesse 
times longer than the resonator cavity length. Thus, from the laboratory frame 
of reference, WGRs also exhibit ultralow threshold SBS lasing [4,16-18] in the 
pW - mW range. This high-finesse resonant enhancement effect is not the focus 
of this paper, but ultimately must not be neglected. 

Here we aim to show that in WGRs the SBS gain experiences additional 
enhancement due to surface optical forces that are ignored by traditional SBS 
theory. To make a fair comparison between waveguides and WGRs we will 
ignore the finesse factor enhancement of the pump and Stokes optical fields. 
Instead, we will only quantify SBS gain as a function of the optical powers for 
pump (Pp) and Stokes (P^) circulating within the WGR. 

In WGRs all three optical and acoustic actors participating in SBS are reso¬ 
nant, which enforces an integer azimuthal mode order. While the modal disper¬ 
sion relationships in waveguides are usually described in uj — k space, in WGRs 
it is preferable to employ n uj — M space description due to this azimuthal dis- 
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Figure 1: Dispersion diagrams showcasing the phase-matching require¬ 
ments for SBS in a WGR. (Top) Stokes scattering in backward SBS, and 
(bottom) in Forward SBS. 


cretization. Here, m or M denotes the integer azimuthal mode order i.e. the 
number of wavelengths around the resonator circumference, and is equivalent 
to the wave momentum k [4,5]. The phase matching requirement for F-SBS is 
then expressed by the equations ujp — oog = and rup — rus = M, presented 
graphically in Fig. 1. In the case of B-SBS, the azimuthal relationship is mod¬ 
ified to Trip + rUg = M. Here, m^, m^, and M correspond to the pump, Stokes 
(scattered) and acoustic modal orders respectively. 

In linear waveguides, the two optical waves can reside in a single waveguide 
mode whose dispersion has been suitably engineered to match the acoustic dis¬ 
persion [28]. In WGRs the phase matching between pump and Stokes signals is 
achieved at avoided crossings of two optical mode families [4], unless anomalous 
dispersion occurs [42]. These optical mode shapes are well studied and easily 
extracted for axisymmetric geometries of WGRs [43,44]. The acoustic modes 
can also be theoretically or computationally evaluated [42] under assumptions 
of circumferential periodicity using a finite element solver [45]. As in linear 
waveguides, the evaluation of SBS gain in WGRs must incorporate the variety 
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Figure 2: Illustration of forces and whispering-gallery optical mode in 
a spherical WGR. For the optical mode, color intensity depicts the magnitude 
of the electric field while the black arrows depict the polarization. Electrostric- 
tion body force appears in the bulk, with color indicating the force magnitude. 
Boundary electrostriction and radiation pressure forces appear at the free sur¬ 
faces, with magnitude and direction indicated by the red arrows. 


of optical and acoustic mode shapes, and their vector nature. 

It is important to consider slight inaccuracies in phase-matching that are 
expected from any real system, which are contigent on the frequencies and 
linewidths of the three modes participating in SBS. The degree of mismatch 
and the resulting SBS gain spectrum can be modeled with the help of detuning 
parameters [20,26] and will certainly affect the achievable gain, but always scales 
with the line-center gain. Thus, in this paper, we consider only the maximum 
line-center SBS gain under the assumption that phase matching is perfectly 
satisfied. 

While not considered here, Kerr and four-wave mixing nonlinearities [46-49] 
are known to appear at very low thresholds in high-Q resonators. These effects 
can become significant in silicon [28] due to its stronger nonlinear response 
than silica, and in other highly nonlinear materials like chalcogenide glasses and 
ferroelectric crystals. 

4 Theory of SBS Gain in WGRs 

Existing full vectorial models of Brillouin gain in linear waveguides [27,30,34- 
36,38] employ the Cartesian coordinate system, which is not suitable for WGR 
analysis. Here we follow the method of ref. [27] to produce a vectorial formu¬ 
lation based on energy and particle conservation for calculating the SBS gain 
in axisymmetric WGRs. The optical force calculations are adapted from the 
theory developed in [30]. Note that geometry differences necessitate notational 
variations with respect to [27]. Our derivation also adds to previous work in [27] 
by directly defining the acoustic amplitude through the optical forces, in order 
to provide a formulation of SBS gain that does not require knowledge of the 
acoustic amplitude. This is possible under the assumption that the acoustic 
and optical force envelope functions have reached steady state. 
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We employ the following conventions: * denotes the complex conjugate, and 
p, 5, subscripts denote variables that correspond to the pump optical wave, 
Stokes optical wave, and acoustic wave, respectively. The cylindrical coordinate 
system is represented by (see Fig. 2) r denoting the radial unit vector, z denoting 
the axial unit vector, and (p is the azimuthal unit vector in the counterclockwise 
direction. The subscripts i, j, /c, and I are used to denote directional variables. 

We begin in a manner similar to [27] by defining the change in Stokes optical 
power Ps with respect to some angular travel 5(j) (units of rad) as follows 

SP 

-j^=GbPpPs-PsPs ( 1 ) 

where Gb (units of rad'^W'^) is defined as the SBS gain per angular travel 
averaged over an acoustic period. Here, Pp is the pump optical power and 
the angular photon loss rate is f3s and has units of rad'^. Note that Pp and 
Ps are defined here as the steady-state circulating optical power within the 
resonator, and are thus related to the input powers via the optical finesse {Pres = 

^oQopt p \ 

The phonon generation rate in this system is related to the change in acoustic 
power per unit angle, which is equivalent to the time averaged overlap between 
the force and the deformation velocity over a differential volume delimited by the 
above Sep. Here, the angle Sep covers an integer number of acoustic wavelengths, 
and the associated resonator volume is 6V. Thus the acoustic power added to 
the system by the force field driving the elastic deformation wave over an integer 
multiple of acoustic periods can be represented by the equation 

SPn= [ (Fn{rp)’MrP))dV. (2) 

Here, FQ{f,t) is defined as the vectorial optical force density (units N/m^) as 
it varies over time and space within the acoustic period-delimited volume, SV. 
Similarly, u^{r^t) is the particle velocity field (units m/s) varying over time 
and space within SV. This particle velocity field is a vector field describing the 
velocity of individual particles in the elastic wave. Additionally, (• • •) represents 
time averaging over an acoustic period, or equivalently, spatial averaging over 
an acoustic period. 

Since particle conservation must apply, the Stokes photon flux generation 
rate (y^ = in the resonator is equal to the phonon flux generation rate 

(70 = ^) • Thus we can equate js = and write 

This expression is the same as that found in [27] with the exception that the 
differential volume is defined by the small angle Sep rather than a small distance 
Sz. Following the same method as [27] we can simplify this expression to 


9b 


ids 

2RQ 


Fv f 

, / Re[P -WWdrdz. 

^/tWsJa 


(4) 
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The new gain variable used in Eq. (4), qb^ represents the Brillouin gain over 
interaction distance (units with the assumption that all interactions oc¬ 

cur near the outermost resonator surface. The force amplitude F{t) and acoustic 
amplitude u(t) have been separated from their complex unitless vectorial mode 
shapes U)^ through the expressions t) = PpPsRe[F{t)F{r^ 
and UQ{fR) = Re[u{t)U{r, The force amplitude is also normal¬ 

ized to remove the dependence on the optical powers. The mode shapes are 
integrated over the half-cross-section of the resonator, A (see Fig. 2), and have 
been defined such that z)\dA = 1 and jj^\U{r^z)\^ dA = 1. These 

are both ‘power-normalized’ integrals in the sense that they are unaffected by 
changes to the optical or acoustic powers. We note here that all normalizations 
above are performed only out to a finite radius, following the standard WGM 
normalizations as used in [50], in order to avoid divergence of the integral. 

We note here that a key difference arises between the gain equation for linear 
systems [27] and the corresponding equation in WGRs (Eq. (4)); a factor of R 
appears in the denominator here due to the conversion from rad'^ to m'^. This 
factor can be thought of as ‘normalizing’ the radial weighting of the field within 
the integral. 

We now remove the acoustic variable ii from Eq. (4) by considering the 
optical forces as a source in the acoustic wave equation. Here we model the 
acoustic wave as a whispering gallery mode which can be written as (r, t) = 
Re[u{t)U{r^ traveling in a circumferential path along the azimuthal 

direction. This acoustic mode, written in terms of the particle displacement field 
flo, is directly related to the particle velocity field u^t introduced previously at 
steady-state hy u = jClu. We can thus write the acoustic wave equation as 

" C Ti ' F 

uq = V • — : + V • — : H-. (5) 

P P P 

Here V- is the tensor divergence, C is the fourth order stiffness tensor, the 
: symbol denotes a tensor inner product. Vs denotes is the symmetric tensor 
gradient, p is the material density, and p is the fourth order viscosity tensor 
[51]. In order to introduce a scalar, phenomenological damping parameter T 
analogous to that used by Boyd [15], we assume the viscosity tensor is directly 
related to the stiffness tensor by = ^C. At steady state the force and acoustic 
envelope F and u are not varying in time, allowing us to arrive at 


uU = 


jRy/RpRs _ 
pQA 


(6) 


We have defined A = T + j(^o ~ as the total loss associated with both 
damping and detuning from the acoustic resonance flo, where T is the phonon 
loss rate in Hz. The line center gain can be calculated by setting the frequency 
Q = CIq. We note that gives the Lorentzian shape of the acoustic resonance as 
expected. Multiplying both sides by hP ^ and integrating over the cross-section, 
we insert Eq. (6) into Eq. (4). We now consider only the line-center SBS gain 



and also introduce the mechanical quality factor, Qo = ^5 to find 

It can be seen that the gain is dependent on the phonon loss rate in the same 
manner as a linear waveguides [15,27]. The nonlinear dependence on the force 
can thus be seen to arise from the relation of the acoustic amplitude to the 
optical force {ii oc F^JPpPg Re[T • hP]dA). 

Finally, the optical force term within Eq. (7) is comprised of electrostriction 
Fes (which has both body and boundary components) and radiation pressure 
Frp, such that F = Fes, body + Fes, boundary + Frp. The line-center SBS gain 
in WGRs is thus given by 


9b = 


2Rp9P 


FES,body / Re[FES,body 

JA 

Fes, boundary j R^[FeS, boundary U*]rdl+ 

Frp j^Re[FRP-W]rdl^ . (8) 


As can be seen, the gain depends on the integral of the phonon generation rate 
over the half-cross-section A of the resonator for the body force, and over the line 
C bounding the half-cross-section for the boundary forces (see Fig. 2 ). Again, 
we note that this gain equation differs from its linear counterpart in [27] since 
the acoustic power is weighted by the radial location R of the wave. 


Evaluating the optical forces 

We now require analytical expressions for the individual surface and boundary 
forces in terms of the optical fields. Here we can assume that the two optical 
signals are confined to a narrow region within a few microns of the resonator 
surface. The azimuthal component of each vector thus points in the direction 
of wave propagation, and the radial and axial components are transverse to 
propagation, similar to a plane wave. 

The electrostriction stress on the material is related to the optical field 
through the photoelastic tensor as shown in [30]. Since silica is an amorphous 
isotropic material, the electrostriction stresses in our cylindrical coordinate sys¬ 
tem can be written through the expression 


CFrr 


2Ep^rFgr^ 
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where the photoelastic tensor P is 
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(10) 


Here, aij denotes the components of the electrostriction stress, n is the effective 
refractive index of the resonator material, the pij terms are the components 
of the photoelastic tensor, the p and s subscripts denote the pump and Stokes 
fields, and the r, z, and <p subscripts denote the directional components. Note 
that aij = aji for the purposes of this paper. 

The electrostriction body force is then obtained from the summed gradients 
of the normal {i = j) and transverse {i ^ j) stresses as shown in [30], which in 
our cylindrical coordinate system are expressed as 


BES, body _ 9 ^ rr 

r — ^rz J ^rcf) 

or oz J. ^ 

(lla) 

pES,bodv___^ -Ter - 7 — er - 

z — o ^ zz o ^ zr J ^ z(b 

oz Or ^ 

(lib) 

BES.body 9 d 

(lie) 

An additional electrostriction boundary force exists at material interfaces, aris¬ 
ing from the difference in electrostriction stresses across the media. Here we 
assume that the resonator is immersed in air, which supports negligible elec¬ 
trostriction stress in comparison to the resonator material (it is zero in a vacuum 
environment). The electrostriction boundary force can then be written as 

pES,boundary ^ 

(12a) 

PES,boundary ^ ^ 

(12b) 

BES,boundary ^ , 

PfP — a(p-pTi-p T a(j)z'biz' 

(12c) 


Here, fir^z are the radial and axial components of the boundary surface normal 
vector h, and are transverse to wave propagation. 

The final contributing force is radiation pressure, which arises from the 
Maxwell Stress Tensor [52]. When considering dielectric systems without free 
charges, as is the case for this study, radiation pressure is localized where the di¬ 
electric constant gradient is nonzero. When the system consists only of domains 
of homogeneous media, radiation pressure exists only on material boundaries. 
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Again considering a resonator immersed in air or vacuum, the boundary radia¬ 
tion pressure force can be written as 

=^-f{e-l)E^^r.El^h (13) 

where e denotes the relative permittivity of the resonator, the n subscript de¬ 
notes the component of the electric field which is normal to the boundary, and 
h denotes the boundary normal vector pointing outwards from the resonator. 
Note that electrostriction and radiation pressure, due to their distinct physical 
origins [53,54], must be included separately as we have done above. 

At first sight the two boundary forces normal to the resonator surface may 
seem mathematically orthogonal to the azimuthal propagation of the acoustic 
wave and unable to contribute to the generation of azimuthal phonons. How¬ 
ever, we remind the reader that these optical forces are a consequence of the 
spatiotemporal interference between the pump and Stokes optical fields. These 
forces oscillate at the pump-Stokes beat frequency, and also spatially propa¬ 
gate with the momentum-vector difference between them. Due to the defined 
SBS phase-matching relationship, this spatiotemporally varying force precisely 
matches the frequency and momentum vector of the acoustic wave, thus co¬ 
propagating with the acoustic wave. Since all surface acoustic waves involve 
some normal displacement of the material surface, boundary-normal forces do 
actuate these waves. In other words, boundary forces have non-zero overlap 
integral with the traveling wave acoustic mode and are able to drive the co¬ 
herent generation of phonons in the system. This interpretation is well known 
from prior work on SBS enhancement in nanoscale waveguides and also in res¬ 
onators [5]. 


Scalar Brillouin Gain 


In order to provide a comparison benchmark for the vectorial SBS gain, we now 
evaluate the scalar SBS gain. Here we assume that a longitudinal acoustic wave 
is driven only by bulk electrostriction as considered in standard SBS theory [15]. 
Additionally, the relation between the optical forces and acoustic wave is further 
simplified by only considering the contribution of the transverse optical fields to 
the electrostriction stress; this simplification reduces the photoelastic tensor to 
a scalar (P 12 ). After making these simplifications, the scalar gain can be written 


as 


M 'n7p\2‘^l Qn V 
rris 2c^ph' ft R 


(14) 


where c is the speed of light and we have introduced p as the spatial over¬ 
lap integral. The spatial overlap integral describes the efficiency with which 
the optical forces drive the acoustic wave and can be evaluated here as 77 = 
{f^Re[T-U*]rdrdz)^. 

For smaller resonators modal confinement is improved, thus the relative over¬ 
lap between optical and acoustic fields tends to increase. As a result, small 
resonators will naturally exhibit a larger SBS gain. However, this trend cannot 
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always be predicted because of the unique shapes of optical and acoustic modes. 
For instance, the symmetric Lamb wave modes considered in the next section 
have sufficient acoustic dispersion variation with decreasing resonator thickness 
to significantly decrease the scalar gain because of the increase in acoustic fre¬ 
quency (gain is quadratically dependent on acoustic frequency Cl as shown in 
Eq. (8)). 

5 Brillouin Gain Calculations and Discussion 

Using the formalism established in Section 4, we now calculate the full vectorial 
Brillouin gain (Eqs. (8)-(13)) and the scalar gain (Eq. (14)) for some cases of 
interest. As we shall see, significant deviations from scalar theory are predicted 
by the vectorial numerical calculations that incorporate the boundary forces. In 
particular, giant SBS gain enhancement up to 10^ greater than scalar theory, 
and significant suppression of SBS gain by as much as 10^ times, are predicted 
in specific cases. 

Cases Considered and Key Assumptions 

Our search here is focused on E-SBS only. This is because B-SBS involves 
extremely high-frequency, short-wavelength acoustic modes that are spectrally 
dense [16,33], preventing accurate identification in experiments on large WGRs. 
However, we stress that the theoretical foundations of Section 4 remain valid 
even for B-SBS. The validity of Eq. (8) for both E-SBS and B-SBS leads to 
the expectation that the trends predicted here (such as scaling of gain with 
resonator geometry) will also appear in the B-SBS case. Overall lower SBS 
gain should be expected for B-SBS due to the inverse-quadratic dependence of 
the gain on acoustic frequency Cl and linear dependence on acoustic Qq, slightly 
counter-balanced by the decreased acoustic mode volume which improves modal 
overlap. Additionally, since E-SBS always has higher gain (phonon lifetimes are 
significantly longer at lower frequencies), we are very interested in how much 
this gain can be enhanced. 

Eour basic WGR geometries (Eig. 3) are considered in this study: solid 
sphere, spherical shell, rectangular disk, and trapezoidal disk. Solid micro¬ 
spheres are easily made on the tips of silica fibers using refiow or arc dis¬ 
charge processes [55]. Spherical shells encompass the microbottle [56] and mi¬ 
crobubble [41] WGR categories. Trapezoidal disks can be fabricated through 
wet-etching [18], while rectangular cross-section disk resonators are fabricated 
through directional ion-etching [57]. 

While many high-order optical WGMs exist in such resonators, here we con¬ 
strain the analysis by considering only four lowest-order TE-like (radially polar¬ 
ized) and TM-like (axially polarized) modes [58,59], all of which are illustrated 
in Eig. 3. The modal subscript numbering convention used here is common to 
WGRs [59] where the first numeral denotes the number of anti-nodes in the ra¬ 
dial direction and the second numeral denotes the number of nodes in the axial 
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Figure 3: Examples of acoustic and optical modes used in this work. 

For acoustic modes, we present equatorial cross-sections and full 3D views of 
the acoustic modes with exaggerated deformation. Red and blue denote the 
magnitude and direction of the displacement. For optical modes, we present az¬ 
imuthal cross-sections of the optical modes in conjunction for each WGR. Four 
lowest optical mode families are examined: TEi^o? TMi^q? TE 2,0 5 TM 2 ,o- 
The color intensity represents the magnitude while the arrows represent the 
polarization of the electric field. TEi^o is radially polarized and concentrated 
in a single antinode. TMi^o is axially polarized and concentrated in a single 
antinode. TE 2 ,o is radially polarized and concentrated in two antinodes, dis¬ 
tributed radially. TM 2 ,o is axially polarized and concentrated in two antinodes, 
distributed radially. 
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direction. We selected the azimuthal mode order rus to place the optical 
frequencies around the 1550 nm telecom band 193 THz) for each case to fa¬ 
cilitate comparison. Pairwise combinations of these simple modes are employed 
as the two (pump and Stokes) optical modes in the computational analysis of 
the SBS gain. 

The third actor in microresonator SBS is the acoustic WGM. Multiple sur¬ 
face acoustic wave mode families are considered in our analysis (Fig. 3), as they 
are known to exist in WGRs [45] and have been experimentally observed as 
well [4,5,32]. Rayleigh waves (radial displacement) are considered in all the 
WGRs that we consider in this analysis, symmetric Lamb waves (extensional) 
are considered in both disk WGR configurations, and anti-symmetric Lamb 
waves (flexural) are considered in the trapezoidal disk WGR configuration. In 
each case, the acoustic mode that satisfies SBS phase-matching between the 
selected optical modes must follow M = rrip — rris as previously described in 
Section 3. However, since each resonator’s size and dispersion determines rup 
near 1550 nm, we risk having no standard point of comparison between the res¬ 
onator geometries. We thus constrain our search to only the lowest transverse 
order acoustic modes, while the azimuthal mode order is fixed at M = 12 as 
a standard comparison point commensurate with past experimental observa¬ 
tions [5,32]. The pump azimuthal mode order rUp is then selected as described 
above, and the availability of a Stokes optical mode with correct rris is assumed. 
This value of M was chosen since numbers in this range have been measured 
experimentally [4,5,32]. 

In the cases where we consider intra-modal scattering (example - the TEi^o/TEi,o 
case) the scattering is taking place between the and {N + azimuthal 

order optical modes having the TEi^o mode shape. Here, N is a large integer 
depending on the circumference of the resonator, while M = 12 is fixed. In each 
calculation, we compute the nearest integer N based on the resonator size. The 
same applies to the inter-modal scattering pairs (example - the TEi^o/TE 2 ,o 
case) where the scattering takes place between the TEi^o mode and the 
{N + TE 2 ,o mode. We also note that the optical forces generated in the 

reversed case of (A/' + M)^^ TEi^o mode and the TE 2 ,o mode are nearly the 
same as long as M ^ A/", since in such a case the mode shapes do not change 
appreciably. This symmetry can be seen in Eqs. (9) and (13) in the evaluation 
of optical stresses and forces. 

In order to set simulation parameters, and to draw comparisons to a com¬ 
monly used optics platform, we have chosen to examine silica WGRs. We note, 
however, that the overall scaling trends observed here are not specific to any 
material platform as they are primarily a geometrical effect. Eor the purpose 
of our numerical computations, the mechanical quality factor is assumed to be 
10^, selected based on experimental observations [4,32]. In smaller resonators, 
the acoustic frequency for an M = 12 mode increases, resulting in a shorter 
phonon lifetime for the same Q-factor. The photoelastic constants for silica [60] 
are set at pn = 0.125, pi 2 = 0.278, and P 44 = —0.073 for evaluating Eq. (9). 

Mode polarization and resonator geometric scale are seen below to be the 
primary drivers for the observed gain enhancement and suppression. Eor in- 
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(a) Total Gain in Spherical WGRs (b) 




—•—Total Gain 
—©— Radiation Pressure 
—B— Electrostriction Body Force 
— 0 — Electrostriction Boundary Force 
-Scalar Gain 


Figure 4: SBS gain computed in sphere WGRs with Rayleigh acoustic 
waves (M = 12) exhibits both enhancement and suppression beyond 
scalar theory, (a) The total gain including all forces, (b)-(e) Gain contri¬ 
butions of individual optical forces if they were considered independently. The 
SBS gain expected from scalar theory is shown for comparison on each plot. 
Note that radiation pressure does not appear in (d)-(e) since it is many orders 
of magnitude smaller than the other optical forces, (f) Legend for sub figures 
(b)-(e). 
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stance, smaller resonators tend to have increased boundary forces and better 
overlap between the acoustic mode and optical forces, resulting in higher gain. 
Brillouin gain is suppressed in cases where radiation pressure force approxi¬ 
mately balances electrostriction. We also observe that the gain predicted by 
scalar theory underestimates the gain predicted by vectorial theory, especially 
for smaller resonator scales where the largest gain enhancements are observed. 

Brillouin Gain in Silica Microspheres 

We first examine the SBS gain in spherical silica WGRs (see Fig. 3) since they 
are extensively used in microresonator studies and are easy to fabricate. The 
simple spherical geometry also helps bring out the relation of the opto-acoustic 
coupling to the resonator scale. Brillouin lasing has been experimentally demon¬ 
strated in microspheres before [4,5,16,20], although the relative contributions of 
radiation pressure and electrostriction have never been analyzed previously in 
the resonator case. In Fig. 4 we present the computed SBS gain in the Rayleigh 
acoustic wave case for sphere radii ranging from 10 — 100 pm using both scalar 
and vectorial calculations. 

We immediately note that while the scalar electrostriction-only theory ap¬ 
proximately predicts the trend shown in Fig. 4, it over-estimates or under¬ 
estimates the full vectorial calculation by an order of magnitude in some cases. 
However, the scalar gain does roughly predict the contributions from each each 
force individually when the optical modes are relatively well confined inside 
the resonator (large sphere diameters), as seen in Fig. 5(b)-(e). Lowered opti¬ 
cal confinement (for smaller sphere diameters) causes the appearance of large 
boundary forces and affects the overlap integral in Eq. (14). These additional 
forces significantly deviate the gain from the scalar estimate since the vectorial 
gain is evaluated through a nonlinear combination of all contributing forces (see 
Eq. (8)). 

As mentioned previously, radiation pressure and electrostriction can oppose 
each other (see Fig. 2) in certain cases. Brillouin gain suppression is seen to occur 
for TE modes (Eig. 4(b)-(c)) when radiation pressure and electrostriction are 
of approximately equivalent magnitude. However, no gain suppression occurs 
for TM modes (Eig. 4(d)-(e)) because the electric field component normal to 
the surface is very small, generating negligible radiation pressure contribution 
as per Eq. (13). 

These results immediately support our main assertion that SBS gain en¬ 
hancement does not require nanoscale or sub-wavelength features as previously 
proposed [27,28,30,31]. SBS gain enhancement can occur whenever optical fields 
are confined along a surface and have good overlap with the acoustic mode. Eur- 
thermore, SBS is not solely driven by bulk electrostriction in resonators as has 
been thought until now, but is instead significantly affected by boundary optical 
forces as well. 
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Shell Thickness (//m) 
1,0 (e) 
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Shell Thickness (fim) 

(f) Legend for subfigures (b) - (e) 

-Total Gain 
-Radiation Pressure 

- Electrostriction Body Force 

- Electrostriction Boundary Force 
-Scalar Gain 


Figure 5: SBS in shell-type WGRs (bubbles, bottles) exhibits very 
significant gain enhancement for the Rayleigh acoustic wave case 
(M = 12). The causes for the gain enhancement are discussed in Fig. 6. (a) The 
total gain including all forces, (b)-(e) Gain contributions of individual optical 
forces if they were considered independently, (f) Legend for subfigures (b)-(e). 
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Figure 6: Gain enhancement in shells with Rayleigh acoustic WGMs 
arises due to boundary and body forces when shell thickness ap¬ 
proaches the optical wavelength. For the body force, red indicates force 
directed radially outwards and blue indicates force directed radially inwards. 
An intra-modal pair and an inter-modal pair are chosen in order to highlight 
the different force distributions in thin (2 pm) and thick (10 pm) shells. The 
force distribution shown external to the resonator does not exist but is shown 
only to indicate the presence of optical fields. 

Brillouin Gain in Silica Microshells 

We now consider gain in spherical shell WGRs of the form shown in Fig. 3. 
Computed gains for the Rayleigh acoustic wave case are presented in Fig. 5 
where the shell thickness is varied between 0.2 — 10 pm using a fixed radius 
of 100 pm. We immediately note giant gain enhancement 10^ x) for thin 
shells, which becomes significant when the shell thickness is on the order of a 
few optical wavelengths (< 5 pm). Since all the modes are confined near the 
outer WGR surface, thicker shells tend to approximate simple spheres, and the 
computed gain also converges to the result for a sphere (compare Fig. 4(a) and 
Fig. 5(a)). A notable feature is the significant gain suppression (^ 10^x) in the 
TEi^o/TE 2 ,o mode pair for the specific case of a 10 pm thick shell. Overall, it 
can be seen that the scalar computation is a poor predictor of Brillouin gain for 
thin shell WGRs. 

For shell WGRs we find that there are two distinct mechanisms of gain en¬ 
hancement. The intra-modal combinations (TEi^o/TEi^o and TMi^q/TMi^q) 
experience gain enhancement due to large optical boundary forces, while the 
inter-modal combinations (TEi^o/TE 2 ,o and TMi^o/TM 2 ,o) experience large 
gains due to the body force. These enhancement mechanisms are evidenced 
in Eig. 6. 
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Enhancement from boundary forces for intra-modal scattering: 

Let us consider the radial component of the combined boundary force for the 
TEi^o/TEi^o mode pair (Eig. 6, top row). When the shell wall thickness de¬ 
creases, the optical modes begin to interact with the inner boundary, and the 
boundary forces thus generated start contributing to the Brillouin gain. Typ¬ 
ically, higher interfacial fields create stronger boundary forces and increased 
Brillouin gain is expected as in the case of nanoscale waveguides. In the mode 
combination studied here, the forces generated on the inner boundary oppose 
the forces on the outer boundary, resulting in a lower contribution to the gain. 

Enhancement from body force for inter-modal scattering: Let us 
now consider the radial component of the electrostriction body force for the 
TEi^o/TE 2 ,o mode pair (Eig. 6, bottom row). Eor this case three regions of 
body force concentration are observed, which alternate between force directed 
radially outwards and radially inwards. These opposing forces approximately 
cancel, and this self-cancellation results in a reduced gain contribution from 
the total body force. Eor thinner resonators, however, the resonator does not 
confine the optical fields sufficiently. The central body force concentration then 
is the primary contributor to electrostriction within the resonator material and 
the above self-cancellation effect is reduced. The gain in thin shells is therefore 
enhanced beyond the value for thicker shells and spheres. The intra-modal 
combinations, on the other hand, exhibit only two opposing force concentrations 
such that the self-cancellation effect persists. 

The significant 10^ x) suppression of gain in the TEi^o/TE 2 ,o mode pair 
for 10 pm thick shells is caused by the comparable and opposite radiation pres¬ 
sure and electrostriction forces in the same manner as described previously in 
spheres. It can be speculated that an exact choice of geometrical parameters 
could also lead to perfect cancellation and thus giant suppression of the SBS 
gain beyond the demonstrated value. 

Brillouin Gain in Silica Trapezoidal Microdisks 

Microdisk resonators with trapezoidal cross-section (Eig. 3) are also widely em¬ 
ployed in photonic microsystems [18,61]. Eor SBS gain calculations, we consider 
three acoustic WGM families: Rayleigh waves, symmetric Lamb waves, and anti¬ 
symmetric Lamb waves (see Eig. 3). We consider disk thickness in the range of 
0.2 — 10 pm using a fixed radius of 100 pm. Anti-symmetric Lamb waves are 
seen to have the largest absolute gain (^ 10^ m'^W'^) while the largest gain 
enhancement (> 10^ x beyond scalar theory) is predicted for symmetric Lamb 
waves. 

Rayleigh acoustic WGMs: These modes are characterized by a large 
radial deflection at the outer edge surface of the resonator. The results of scalar 
and full numerical computation in this case are presented in Eig. 7. Unlike 
solid spheres and spherical shells, we predict that the gain of mode pairs of 
both TE and TM polarizations is primarily driven by radiation pressure. The 
axially polarized (TM) modes are able to drive a radially directed acoustic wave 
by means of radiation pressure because the angled edge of the trapezoidal disk 
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(f) Legend for subfigures (b) - (e) 
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Figure 7: SBS in trapezoidal disk WGRs with Rayleigh waves (M = 12) 
is well approximated by the scalar gain. TE mode pairs exhibit a sudden 
rise in total gain in disks of thickness 2 — 5 pm. (a) The total gain including 
all forces, (b)-(e) Gain contributions of individual optical forces if they were 
considered independently, (f) Legend for subfigures (b)-(e). See also Fig. 8. 
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Figure 8: The sudden rise of gain in trapezoidal disk WGRs is at¬ 
tributable to the location and polarization of the optical modes. The 

lowest order mode of each polarization is shown in thin (2 pm) and thick (10 
pm) disks. The red arrows indicate the direction and magnitude of the in-plane 
electric field. 


breaks symmetry. 

For thick disks TM polarized optical mode pairs show much larger radiation 
pressure contribution than TE polarized mode pairs. This is because the TE 
mode boundary-normal E-fields at the outer edge are very small (Eq. (13)), as 
is simulated in Eig. 8 (bottom left). In contrast, the TM modes exhibit larger 
boundary normal fields (Eig. 8 bottom right). Eor thinner disks (< 2 pm) that 
are comparable to the optical wavelength, the TE optical modes have a larger 
boundary-normal E-field component at the outer edge, as is simulated in Eig. 8 
(top left). This increases the radiation pressure contribution for the smaller size 
scales, with the Brillouin gain transition for Rayleigh modes being easily visible 
in Eig. 7(b)(c). 

Symmetric Lamb wave WGMs: These modes are characterized by mir¬ 
rored axial deflection along the thickness axis of the disk [51] as illustrated in 
Eig. 3. The scalar and full numerical computation results for this case are shown 
in Eig. 9. Eor very thin disks the distinction between symmetric Lamb waves 
and longitudinal waves is no longer apparent [62], so we only calculated gains 
for thickness ranging over 2 — 10 pm. 

Symmetric Lamb modes exhibit periodic compression and rarefaction, much 
like a longitudinal wave in a bulk material, and intuitively ought to exhibit 
high SBS gain. However, the regions of axial compression in these acoustic 
WGMs are accompanied by radial expansion of comparable magnitude. As a 
result, self-cancellation can occur in the overlap integral of the compressive-only 
electrostriction force with this mode shape. Similarly, the boundary forces drive 
the acoustic mode axially via the upper and lower faces, but oppose the mode 


21 






















































Total Gain in Trapezoidal Disk WGRs with Symmetric Lamb Waves 
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(f) Legend for subfigures (b) - (e) 
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Figure 9: SBS in trapezoidal disk WGRs with symmetric Lamb waves 
(M = 12) has low total gain but shows enhancement over scalar the¬ 
ory. Disks of thickness smaller than 2 [jim were not tested since symmetric 
Lamb waves transition to longitudinal waves in thin plates. Acoustic dispersion 
of these symmetric Lamb modes is responsible the downward trend of scalar 
gain for thinner disks, (a) The total gain including all forces, (b)-(e) Gain 
contributions of individual optical forces if they were considered independently, 
(f) Legend for subfigures (b)-(e). 
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(f) Legend for subfigures (b) - (e) 
-Total Gain 
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Figure 10: SBS in trapezoidal disk WGRs with anti-symmetric Lamb 
waves (M = 12) exhibits extremely large gain, and some enhancement 
over scalar theory. Significant reduction of the electrostriction body force is 
observed for TE modes in the 1 — 5 [im thickness range due to a transition in 
the optical mode shape (see Fig. 8 and main text), (a) The total gain including 
all forces, (b)-(e) Gain contributions of individual optical forces if they were 
considered independently, (f) Legend for subfigures (b)-(e). 
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Figure 11: SBS in rectangular disk WGRs with Rayleigh waves (M = 
12) exhibits significant gain suppression in TE modes due to cancel¬ 
lation of radiation pressure and electrostriction forces, (a) The total 
gain including all forces, (b)-(e) Gain contributions of individual optical forces 
if they were considered independently, (f) Legend for subfigures (b)-(e). 
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radially along the outer edge. Rayleigh and anti-symmetric Lamb waves are 
not affected in this way as they primarily involve displacement along a single 
direction. Thus, Brillouin gain in disk WGRs with symmetric Lamb waves is 
not significantly enhanced by thinning the disk. 

Another notable trend observable in Fig. 9 is the decreasing scalar gain as 
a function of decreasing disk thickness. Symmetric Lamb waves have higher 
acoustic frequency on thinner substrates, which in this case should decrease the 
gain quadratically according to Eqs. (8) and (14). 

Anti-symmetric Lamb wave WGMs: These modes are characterized 
by simultaneous upwards or downwards defiections of the upper and lower res¬ 
onator surfaces [51]. The results of scalar and full-vectorial computation are 
presented in Fig. 10. We note that very large SBS gain is attainable in thin 
disks, and is relatively well predicted by scalar theory. Here, the gain increase 
primarily occurs due to the decreasing acoustic frequency of the anti-symmetric 
Lamb wave; SBS gain in Eq. (8) increases quadratically as acoustic frequency Q 
decreases. This trend matches the known acoustic dispersion for anti-symmetric 
Lamb waves on thin substrates [51]. 

Additionally, we observe significant reduction of the electrostriction body 
force contributions for thin disks, in the regime where the thickness is compara¬ 
ble to the optical wavelength (Eig. 10(b)-(e)). This effect is caused by the disk 
edge, where the optical electric field orientation transitions from tip-oriented 
(> 5 pm) to disk-oriented (< 1 pm) as simulated in Eig. 8. In this transition 
region the optical modes change shape such that the electrostriction body force 
is not able to drive the acoustic mode efficiently. 

Brillouin Gain in Silica Rectangular Cross-Section Microdisks 

We now examine Brillouin gain in microdisk resonators of rectangular cross- 
section, as illustrated in Eig. 3. In this case we only analyze Rayleigh and 
symmetric Lamb-wave acoustic modes. Here, the disk thickness is varied be¬ 
tween 0.2 — 10 pm using a fixed radius of 100 pm. Overall, we see that the gain 
predictions for rectangular disks are similar to their trapezoidal counterparts 
with some deviations induced by the absence of the angled outer edge. 

Rayleigh wave WGMs: The results of scalar and full-vectorial computa¬ 
tion are presented in Eig. 11. We note that the gain for axially polarized (TM) 
mode pairs is similar to the computed values from the trapezoidal disks, i.e. lies 
in the 10“^ — 10“^ m'^W'^ range. In contrast to trapezoidal disks, however, 
here the radiation pressure contribution to the TM cases is negligible due to 
absence of the angled outer edge. 

Radially polarized mode (TE) cases do show much higher radiation pressure, 
on par with electrostriction forces, however the total gain is seen to be much 
lower than the TM case. This is because in the TE case these forces oppose each 
other resulting in cancellation of contributions. As a result of this, significant 
gain suppression(lO^x) is expected for the radially polarized inter-modal pair 
(TEi,o/TE2,o). 
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Total Gain in Rectanguiar Disk WGRs with Symmetric Lamb Waves 
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Figure 12: SBS in rectangular disk WGRs with symmetric Lamb waves 
(M = 12) exhibits low total gain. Significant suppression of SBS with the 
higher order TE mode pair is also noted. Acoustic dispersion of these symmetric 
Lamb modes is responsible the downward trend of scalar gain for thinner disks. 

(a) The total gain including all forces, (b)-(e) Gain contributions of individual 
optical forces if they were considered independently, (f) Legend for sub figures 

(b) -(e). 
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Symmetric Lamb wave WGMs: The results of scalar and full-vectorial 
computation are presented in Fig. 12. We predict the gain to be similar to that 
in trapezoidal disk WGRs. Again, the gain decreases in thinner disks due to 
the acoustic dispersion. 

Anti-Symmetric Lamb wave WGMs: In a thin rectangular cross-section 
disk, the optical modes are mirror-symmetric relative to the plane cutting the 
middle of the disk. The overlap integral (Eq. (8)) of the resulting optical forces 
with anti-symmetric Lamb wave modes yields zero. Thus, we do not analyze 
this situation in this study as the computational results appear simply as noise 
of numerical calculation. In stark contrast, as seen above, the Brillouin gain is 
extremely high in trapezoidal cross-section disks since this mirror-symmetry is 
broken. 


6 Conclusions 

In this study a full vectorial model for computing the SBS gain associated 
with surface confined acoustic and optical modes in WGRs is developed. The 
model is applied to a variety of common WGRs of spherical, shell, trapezoidal 
disk, and rectangular disk geometries. The overall SBS gain is evaluated by 
analysis of individual contributions from bulk (electrostriction) and boundary 
(electrostriction and radiation pressure) optical forces. Using this model, we 
support our hypothesis that giant enhancements of SBS gain generally take 
place where optical and acoustic fields are simultaneously confined at a free 
surface or material interface, and do not require nanoscale features. 

Our computational results show that significant boundary forces are present 
in microscale WGRs, causing both amplification and suppression of the SBS gain 
of the system. This is contrasted against linear waveguides where nanoscale di¬ 
mensions [27,28,30,31] or slotted geometries [38] are essential to achieve similar 
effects. As expected from previous studies, decreasing the resonator dimensions 
generally increases the boundary forces resulting in overall gain enhancement. 
We further note that specific WGR modal configurations (mode shape and po¬ 
larization) and geometries can be used to enhance or suppress the SBS gain by 
very large factors, even without considering resonance effects. In particular, we 
predict extremely high gain enhancement (by a factor of 10^) in shell WGRs, 
and gain suppression (by a factor of 10^) in rectangular disk WGRs, for partic¬ 
ular sets of modes. It is notable that extreme SBS gain up to 10^ m'^W'^ can 
appear in trapezoidal disks without relying on any resonant enhancement effect. 
We expect these effects to be ubiquitous across WGRs of any dielectric since the 
enhancement due to boundary forces is a geometric effect. The specific ratio 
of contribution from boundary forces and bulk electrostriction, however, will 
be modified by parameters of the particular material, including permittivity, 
photoelastic coefficients, and acoustic parameters. 

Finally, we remind the reader that a further amplification factor of optical- 
finesse-squared, once for the pump power enhancement and once for the Stokes 
path length enhancement, applies to all the gains that we compute in this study 
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(see Section 3). This factor applies only to resonant cavities (including linear 
resonators) and can approach 10^ x 10^ in ultra-high-Q microresonators. Ex¬ 
tremely high SBS gain of 4 x 10^ was recently experimentally measured 

in a microsphere WGR [20] having 10^ finesse factor for both optical modes, 
suggesting a ‘bare’ gain of 4 x 10“^ m'^W'^. This number is consistent with 
the enhanced gain (Fig. 4) presented in this work and is one order-of-magnitude 
higher than the prediction of scalar theory. Extreme Brillouin gains approach¬ 
ing 10^^ m'^W'^ may thus be reachable in ultra-high-Q trapezoidal microdisk 
WGRs (Fig. 10), which is 10^ times greater than the highest predicted [27,38] 
and observed [29] gains in linear waveguides. 
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